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PROBABILISTIC REPRESENTATION OF FUNDAMENTAL 
SOLUTIONS TO || = K™f^. 



ENZO ORSINGHER AND MIRKO D'OVIDIO 



f**^ Abstract. For the fundamental solutions of heat-type equations of order n we 

f\l give a general stochastic representation in terms of damped oscillations with 

J generalized gamma distributed parameters. By composing the pseudo-process 

jrt ^n related to the higher-order heat-type equation with positively skewed stable 

^H r.v.'s Tj , J = 1, 2, . . . , n we obtain genuine r.v.'s whose explicit distribution 

^^ is given for n = 3 in terms of Cauchy asymmetric laws. We also prove that 

■^ Xz{T\ (. . . {T'Kt)) . . .)) has a stable asymmetric law. 



Ph 1. Introduction 

(-H 

'^j The problem of studying the form of fundamental solutions of higher-order heat 

J5 equations of the form 

' ' (1.1) "^{x,t)^K„,'^^{x,t), xeM, i>0, m>2 

' with 

^^ where 



±1 if m is odd 

rn . . " 

f~^ has been tackled in some particular cases by mathematicians of the caliber of Bern- 

p\| stein [3]; Levy [3; Polya [TT] and Burwell |1|. By applying the steepest descent 

T— I method some recent papers by Li and Wong [5' , Accetta and Orsingher T , Lachal 



[B] have explored the form of the fundamental solutions of equation ( 1.1 ). The aim 



• '-j of this note is to give an explicit representation of the solutions to ( 1.1 ) for the case 

rS where the order of the equation is odd, alternative to the inverse Fourier transform 

"^ (1.2) u,,{x,t) ^^ f e-^-+*(-/')'"-™d/? 

and capable of representing the sign- varying behavior of the fundamental solutions 



to ( 1.1 1. Our result is that the fundamental solutions to ( 1.1 1 have the probabilistic 



representation 



(1.3) U2n+i{x, t) = —E L-b'^-G'-^'i^/^) sin (a„a;G2"+i(l/t)) 



XTT 



X eu, t > 
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in the odd-order case, and 



(1.4) 



U2n{x,t) 



1 



-E 



sm 



{xG'"{l/t)) 



X e 



t > 



for the even-order case. In (1.3 1 and ( |1.4[ ) by G'^{t) we denote the generaUzed 
gamma r.v. with density 



g"<{x,t) =7- 



t 



exp 



t 



x,t > 0, 7 > 0. 



The parameters an,b„ appearing in (1.31 and (1.4 1 are 



cos 



2(2n+l) 



bn = sin 



2{2n + l) 



Results (1.3 1 and (1.4 1 show that the fundamental solutions have an oscillating 
behavior which has been explored in several papers by many researchers. In our 
view our result represents a concluding picture of the solutions to higher-order heat 
equations. For all values of the degree n of the equation (1.1) we have solutions 
which have the behavior of damped oscillations where the probabilistic ingredients 
(the generalized gamma or Weibull-type distributions) depend only on n € N. An 
alternative universal representation of the fundamental solution in the odd-order 
case reads 

(1.5) U2n+iix,t) = ^(^^^1 -sin^^ni 



k=l 



"+^tj k\ 2n + l 



2n + l 



Functions U2n+i display ascillations which fade off as the degree 2n + 1 of the 
equation increases. A special attention has been devoted to third-order equations 
where we have that 



(1.6) 
(1.7) 



U3{x,t) 



1 



Ai 



\/3t 

3t r°° 



X 

2 sm 



X e 
'V3 



t>0 



xy y^e *^ dy 



^ OO 

-E 

TTX ^ — ' 
fc=l 



X ' 



1 . T^fc 



In the fourth-order case (biquadratic heat-equation) in Orsingher and D'Ovidio [10] 
we have shown that 

U4{x,t) 



1 

'2n 



_y21 

e 22 cos (xy) dy 



1 



27rV2tV^fr^ (2fc)! 



1/2 ^ 



V2\x\ 

" il/4 



2fc 



r 



In a recent paper we have shown that the composition of an odd-order pseudo- 
process X2„+i with a positivly skewed stable r.v. T i of order 2n+i yields a 
genuine r.v. with asymmetric Cauchy distribution, that is 

tcos -. 
(1.9) Pr{X^ (T^ {t)) e dx} = —. ^-^^^ .dx. 





'•^""2(2n+l) 


TT 


[(x + tsm2(2:+i))2-ft2cos%(2^^,)] 
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For n = 1 from (1.9) we can extract a very interesting relationship for the Airy 
function which reads 

Pr{Xi{Ti{t)) e dx}/dx^ 



1 



-.Ai 



t 1 
'"is J s ^ 



Ai 



ds = 



V5 



>o Vis \^/3sJs^3s \</3sJ 2Trx^+xt + t^' 

We show here that the TTi-times iterated pseudo- process (with rf_j_ , j = 1,2, 

2n + l 

independent stable r.v.'s) 

z,n{t) = X2n+i(ri^(Tl^(. . . (r!!^(t)) . . .))), t > 



2ti + 1 2n + l 



is a stable r.v. of order 
(1.10) 



(2n+l)'^ 



rr with characteristic function 



-i|/3|( 



cos ■ 






2(2n + l)™ ■ ^|/3| 2(2n+l) 
We have also explored the connection between solutions of fractional equations 

d^u du 



0, x>0, t>0 



with the solutions of higher-order heat-type equations (1.1 ) for a 

2. Pseudo-processes 



-, m £ 



Some basic facts about the fundamental solutions of higher-order heat equations 
had been established many years ago essentially by applying the steepest descent 
method. In particular, Li and Wong 8\ have shown that the number of zeros is 
infinite for solutions to even-order equations. The steepest descent method was 
applied by Accetta and Orsingher [T for the analysis of the third-order equation. 
The oscillating behavior of the solutions of higher-order heat-type equations is 
confirmed by our analysis. Furthermore, for the odd-order case our results show 
that the asymmetry of solutions decreases as the order 2n + l increases. The result 
of Theorem |2.1| below shows that solutions of all odd-order heat equations can 
be constructed by means of damped oscillating functions with gamma distributed 
parameters. 

We pass now to our principal result. 



Theorem 2.1. The solution to 

du 

(2.1) 



u{x, 0) = S{x) 



X e 



t> 



is given by 

(2.2) U2n+lix,t) 



-E 



-b„xG^ 



'(!/*) Sin (a„xG2"+i(l/i)) 



where, Vi > 0, the r.v. C^lt) has the generalized gamma distribution 



g^{x,t) = 7- 



r7-l 



exp 



x,t > 0, 7 > 0. 



and 



2(2n+l)' 



2(2n+l) 
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Proof. We start by evaluating the Fourier transform of (2.2) 
(2.3) f e'^''u2n+i{x,t)dx 

3'^-dx r (2" + ^)t ^-b„.w sin{a^xw)w^''e-*^"''^'dw 



= (2n + l)i/ w-'^e- 



,,2^+1 / e 



^if^x+ianXW—bjiXW A0x — ianXW — bnXW 



2'Kix 



-dxdw 



= {2n + l)i / «;^"e-*'" [H0{w{a„ - ib^)) ~ Hp{~w{a^ + i&„))] dw 
Jo 

where in the last step we used the integral representation of the Heaviside function 



Hyix) 



1 

2^ 



giyw 



-dw 



iw 



1 

2^ 



-lyw 



-dw 



By a change of variable, the Fourier transform (2.3) takes the form 

(2n + l)i 



e'-P''u2n+iix,t)dx 



(a„ - ^6„)2«+l 7o 

(2n + l)t 
~ (a„ + z&„)2»+i Jo 



,,2"p-*(j^^^) 



Hp{w)dw 



w2"e-H^r;rTWJ Hp{~w)dw 



(2.4) 



=i(2n + l)t/ w;2"e-"'"'"^'H^(w)dw 



+ i(2n+l)t/ w^"e*"" Hp{-w)dw 
Jo 

/ + 00 
u;2"e~'*"'"'^'iJ^(u;)dw 
-OO 

--i{2n + l)t u;2"e-'*"''"^'dw. 

J/3 



In the above steps we used the fact that 

(a„ + i6„)2"+i = z, and (a„ - ifo„)2"+i 



The integral (2.4) can be performed in two different ways. First we can take the 

Ht{2n+ \)dtdw 



Laplace transform 

g-pt ( j e''^^u2n+i{x,t)dx]dt 



P 



^2n I g_(;,+.™2„+ljj^. 



i{2n + l)w2"dw 
p (/x + iu;2"+i)2 

1 



/■OO 

Jo 



This shows that 



+ 00 



e'P''u2n+i{x,t)dx = e" 



-it/3'' 
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We can arrive at the some result by means of the following trick 

i{2n+l)t I u;2"e 



e U2n+i(a;,t)(ia; = lim 
lim 



2n„-itM)^'*+^-/ito^"+^ 7 



it 



JfS 



Ai->0 11 + it 



We have thus shown that the Fourier transform of (2.2) coincides with the Fourier 
transform of the solution to the Cauchy problem (2.1). D 



For the special case of the third-order heat equation we have the following result. 
Theorem 2.2. The solution of the Cauchy problem 

xeR,t>0 



(2.5) 

can be written as 

(2.6) 

(2.7) 

(2.8) 



du _ _a_u 

dt ~ dx^ ' 

u{x,0) — S{x) 



U3{x,t) = E 

XTT 







e 2G^(i/*)sinf^G3(lA) 



e 2 sm —xy I y^e '^ dy 



-rj=Al —= 



xeR,t>0. 



Proof. It is convenient to work with the following series expansion of the Airy 
function (see Orsingher and Beghin [5', formula (4.10)]) 
3-2/3 



(2.9) 



Aiiw) 



^(3V3^) 



i/3„,AfcSin(2^(fc + l)/3)^/fc + l 



k>0 



k\ 



3-2/3 



IT ^^^ ""' {k-l)\ ^ 

k>l ^ ' 

^^^ /3 ,sm(27rfc/3)^/fc^^ 



If we expand the function 



(2.10) 



^{x,(j)) =e^™"'^sin(j;sin(/)) ^ e^''°^'t'- 



2i 



S2, j-k ^i(t>k _ g- 



E 

fe=0 



1 
" Yi 

-i4>k 



2i 

■^ {xe^^f _ ^ (xe 



.k=0 



fc=0 



k\ 



00 u 
X 



k\ 



2i 



E 

fc=0 



fc! 



we establish a relationship which is useful in transforming (2.9) as 
Miw) =^ y(3i/3^). Bin(2^fc/3) ^ /fc ^ ^^ 



k>l 
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^TT^' fc! io 



fe>i 



fe>i 

WTTJo 



g-^g-(32)i/^W2 



3'^' 1/3 \ ^ 

sm I ——z ' w ] az. 



Now wc write 
1 



(2.11) 



3t 



(3i)i/3^HM^j =^X -^exp^-z^i-:^jsin^— zxUz. 



/ 3 za;x . /3I/2 \ 



From (2.2), for ?i = 1 (7 = 3), we write 
3t '■°° 



U3(x,i) 



TTX 



in I xz cos — j 



6 sin I xz cos — z e dz 



1 



= (z^i = y) = — / e-ye 



sm a; 4/ - du. 

2 V U ^ 



This proves, in a different way, that 



"^^"''^^^^H^ 



D 



Theorem 2.3. We can write the fundamental solution U2n+i in the following al- 
ternative form 



1 



(2.12) W2n+l(x,t)= V 

irnr ^ — ^ 



a; \ 1 niik 

sm r 1 



■nx^ \ '■'^^tj k\ 2n + l 



2n + l 



Proof. From (2.2) we have that 

(2n+l)^ p -^ysin^,^^^ . ( ^ ... 

W2«+i(a;,i)= / e 2(2„+i)sm xycos-— — -—^\V e 

TTX Jq \ 2(2n+ 1' 



2".-*^'"^'dy 



(2n+l)i 

TTX .70 



e " (2n+i) gjjj 2;ysm— lye" ay 



(2n+l) 



=(by(2.10)) 



^a: io ^n fc' 2n + l ^ 



—X \ ^ ■ ( nuk 
TTX^-'\^^iJ fc! ^™ V2?i + i; V ' 2n+l 



^ 00 

-E 



D 



Remark 2.4. We note that 

W2n+l(0,t) 



1 f rnr \ f 1 

sin F 1 



n— >-oo 

2n + 1 J TT 



TrissTT \2n + l 
Remark 2.5. We are able to evaluate the semi-infinite integral 

U2n+l{x,t)dx 
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by means of the representation (2.2). Indeed, by considering that 
(2.13) 



sm (Ax) ax = arctan — 

X ^ ' \B 





vr [A 

= --arccot(^- 

(see Gradshteyn and Ryzhik [5, formula 3.941]) we get that 



U2n+l{x,t)dx = / 

Jo XTT 

TT 



.2 

TT 

2 
1- 



E 
arccot 

arccot 
1 



-6a;G^"+"(l/t) 



sin ('aa;G'2"+i 



(lA)) 



dx 



(i) 



TT 1 



2 2n+l 



2 2n+l 



(2.14) 

^ ^ "> 2n+l 

which is in accord with Lachal [6l formula 11] 

Theorem 2.6. The solution to 
du 

with initial condition u(x, 0) — 5{x) can he written as 

1 



(2.15) 



«2n 

(-1)"+^^^, xGM,i>0 



(2.16) 



U2n{x,t) = -^E 
■KX 



3in(a;G2"(l/t)) 



Proof. The solution to (2.15) is given by 

U2n[X,t) = - 
TT 

By integrating by parts we get that 



,-t)9^ 



, , sin /3a; ,02 



TTX 



+ 



2nt 



cos f3xdf3. 



32n~l-tp^ 



sin /3x dj3 



and this concludes the proof. 



Remark 2.7. The solution (2.16) takes the following values in the origin 

1 



u2„(o,f) = ^ii;G2"(iA)= ^^^^^^ 



r 1 



t \ n— >-oo t 

2n/ ^ tt' 



D 



Clearly, in force of the symmetry and by considering formula (2.13), we obtain that 

1 

U2n{x,t)dx = -. 
^ 

Remark 2.8. It is well-known that the solution to the fractional diffusion equation 

d^ u \2 d u 



(2.17) 
with 



g^-A2fJ, xeM, f >0,j/e (0,2] 
u(x, 0) = 5{x) 



Mt(a;,0)=0 for i/e(l,2] 
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is given by 



2AW2 -2'i-2 



Af-^/s 



^ivl2 Z^ fc! I \tvl2 J 



2AW2 ^ fc! V AW2 ; r (1 - f (1 + fc)) 



^f^(-J^)r(i(i.^-))..»(T<'-^') 



27rAW2^^/c! V AW2 

i:=0 



The folded solution to the equation (|2.17|) reads 

u^{x,t) = 



2ui,{x,t), a; > 
0, a; < 



and for ly = 2a, a € (0, 1), A = 1 becomes 
(2.18) ..(a.,0=^E(-^) 



X \'- Tla(r + 1)) 



This represents a probability density of a r.v. X{t) on the half-line (0, oo) which 
can be expressed in terms of positively skewed stable densities. 

/ X t \ /i"+^ \ ai"+^ 

Pr{X{t) e dx} =q^{x,t)dx = (^- = — j = g„ (^^^j 't^'^'" 



^Em(-^J rKl + fc))sm(^a(l + fc))^^^dy 



7rt' 
at 



id ■^-^ 1 



fc=0 



fc! V 2/" 



r(Q;(l + k)) sin (7ra(l + k)) dy 



CX3 ^ 

E -n^H^ + k))y-''^''+^'>-H^+^ sin (7ra(l + k)) dy 



a 



TT ^ — ' fc! 

fc=0 



=Pa{y,a,t)dy 
where in the last step the expression of the stable density 

oo _, 

(2.19) Pa(x,a,l)^-Y'—r(a(l + k))x-°'^^+'''>-^sin(na{l + k)), xeR+. 

TT ^-^ k\ 
k=0 



The calculations above show that the r.v. X{t) with distribution (2.181 can be 
expressed as 

where Y{t) is a positively skewed stable-distributed r.v. of order a G (0,1). In 
other words the stable law of Y(t) is related to the folded solution of the fractional 
diffusion equation X{t) in the sense that 



Y{t) =■ i ' ^ 



X{t) 



1/q 
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This is because 






{x,a,t)dx. 



We give also the Laplace transforms with respect to time t and space x of qa{x,t) 



(2.20) 



Jo 



e ^qa{x,t)dx = — 



1 / 1 



where in the last step, formula 

1 



£^aa(-Ar), t>0 



X \x J 

has been applied (see formula (5.1) of Beghin and Orsingher [2]). Furthermore, 

/>oo 

(2.21) / e-'^*g„(x,i)dt = /i"-ie~^^°, a; > 0. 



Formulas above help to check that qa satisfies the fractional equation 

/„ „„x d°"u du 

2.22 + = 

^ ' c»i" dx 

with iniital condition u{x, 0) — 5{x) where -^ is the Caputo fractional derivative. 
Remark 2.9. The double Laplace transform of ( 2.22[ ) reads 
(2.23) /i"/:(A,M)-/i"-i=/:(A,/i) 

where 



/o 



-^t 



-Xx 



qa{x,t) dxdt. 



Since 

(2.24) 



C{\,y) 



M 



/i" + A 



the i-inverse Laplace transform is (2.20) while the a;-inverse Laplace transform can 

> 



be obtained from (2.24|) as 

(2.25) 



Q-l 



Ai" + A 



^^ 



Remark 2.10. For a = 1/2, the formula (2.19) yields 



Pi/2(x;l/2,l)=— ^(-1 



,r((r + l)/2)_(::+i,_i 



X 5 -^sin(7r(r + l)/2) 



27r 



—^, — -^.T '^ 5 sm Trfc + - 
(2fc)! V 2 



E 

fc=0 

1 - r(fc + i) _, . / ^ 



r(fc) 

(2/c-l)! 



X sin (Trfc) 



27ra;3/2 ^ (2fc)! 

. 1 f ,_i,>£(i±i).-. 



"27ra;3/2 



fc=0 



(2fc) 



10 



ENZO ORSINGHER AND MIRKO D'OVIDIO 



'27ra;3/2 ^ (2fc)l 



r(fc) 



^ fc=0 



E 



( 1) 2^^''x~'' 



(2.26) 



V47rx3 V 4a; 



The result (2.26) shows that the stable law (2.19) for a. — 1/2 coincides with the 
first-passage time of a standard Brownian motion through level l/v2- 



Remark 2.11. Two solutions to the third-order p.d.e 

du d^u 

dt dx^ 



(2.27) 

are given by 



p{x,t) 



1 .Y 2; 

-=Ai —= 



and 



Indeed, we have that 



and 



,{x,t)='^^^M 



X 



^-p{x,t). 



dq 
di 



X X UJ) 

{x, t)=- -^p{x, t) + - — {x,t) 



dq , , 1 , , X dp , , 

-ix,t)^-p{x,t) + --{x,t), 



dx"^ 
d^q 



t dx 
3d^p, 



tdx^ 



d^p 






tdx^ 



t dt 



and therefore 



dq d'^q 1 / d^p 



By observing that 
dx"^ 



1 f X 
Wt^VWt 



- •«""<" 



3i^ 



'it 



and the fact that 

y"{z) - zy{z) = 
that is, y satisfies the Airy equation, we get that 

dt^^'^dx^^'^ 
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By recursive arguments it can be shown that 

'x\™- 1 



(2.28) 



fm{x,t) = i^-j 



it \^/it) 



is also a solution to (2.27) for m > 1. 



Remark 2.12. We have shown in a previous paper that the r.v. 

(2.29) Z(t)=X3(ri(t)) 

obtained by composing the third-order pseudo-process X^ with the stable subordi- 
nator Ti with distribution 



(2.30) 

possesses Cauchy distribution 



Pr{Ti (t) e ds] = ds-^=Ai 
■' s -^ 



t 



/3s 



s,i>0 



Pr{Z{t) e dx} =dx r -^M (^) \^M (^) ds 

V3t 



(2.31) 



=dx 

—dx 



2tt{x^ +xt + f^) 

V^ t 

2^ (a; +1)2 + 1^2 



Result ( 2.31 ) shows that ( 2.29 ) is a genuine r.v. .The characteristic function of (2.31 ) 
is clearly 

/-l-oo 
e^P^Pr{Z{t) e dx} 
-OO 



g-^t|/3|-»|/3^ 



We have now the following generalization of the previous result for the compo- 
sition of the pseudo-process X^ with successively composed subordinators of order 



Theorem 2.13. The r.v. 

(2.33) Z„(t)=X3(Ti(...(rr(i) )...)) 



t\p\ 3"-i (^cos ^-^ + i sgn{!3) sin ^-^) 



with T{, j — 1,2, ... ,n, independent, positively skewed r.v.'s with law (2.30) has 

3 ' ' 

characteristic function 

Ee^PZ^i') =exp 

(2.34) =exp 

Proof. We first observe that 
Pr{Zn{t) e dx}/dx 



-t. 



cos 



2-3 



vr / 

— [1 + i sgn{P) tan 



2-3'^ 











/3si 



-.Ai 



'' ' -.m' '' 



3si/ si ^351 



/Ssi 



t 1 



' OSfi 



-.Ai 



/ oSf^ 



Yldsj 



has Fourier transform 
e 



'^^Pr{Zn(t) e dx} 
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\-rP^s,-!2^^J S2 



Si ^^3^ 



/3si 



Sn \ OSfi 



-.Ai 



y]i{d^. 



/3s: 



j=i 



S2^3i^ 



3S2 



< 1 



.Ai 






i=2 



= exp -t(e*Ts5n(/3)|/3p)'" 

1 - 

= exp -t(e*t''9"('^)|/3|3V" 

= cxp -~t\P\ 1^ (^cos ^-^ + i sgn{(3) sin ^^) 
We observe that the characteristic function of a stable r.v. Sa{t) can be written as 
(2.35) £;e*/3S„(t) 

= exp 



^g-ti/jre-^^ 



where 



-ai|/3r(l-z0^tan^ 



tan ^ ^ 



and cr = cos ^ > 0. In our case a ~ „„_i , j/ = tt;^ and therefore cr = 



cos 



2-3" 



and 



(2.36) 



tan^ 
tan;^^ 



and since 6 ^ ±1, the r.v. Zn is spread on the whole line with parameter of 
asymmetry equal to (2.36). D 



Remark 2.14. We note that by adjusting the derivation of (2.34) we can obtain 
result (1.10 1 of the introduction. 

Remark 2.15. The positively skewed stable r.v. Ta{t), t > 0, a e (0,1), with 
Laplace transform 

has characteristic function 

_ggi/3T„(t) ^£;g-(-«/3)T„(t) 



cxp 


-i|/3|"e-*'^"S"('') 






exp 


-t|/3r(cos™-zs5n(/?)sin™)" 


exp 


-il/^rcos^fl- 
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and therefore with asymmetric parameter 9 = +1 and a = cos ^. 
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Figure 1. Profiles of odd-order solutions, m = 2n + 1 
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